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We report on a first-principles investigation of the electronic structure and of the magnetic prop- 
erties of the quasi- two- dimensional Mott insulator SrCu2 (603)2- Based on the hopping integrals 
and Coulomb interactions calculated with LDA and LSDA+U, we provide a microscopic explana- 
tion of the symmetric Heisenberg and antisymmetric Dzyaloshinskii-Moriya exchange integrals of 
SrCu2 (603)2. The intra-atomic exchange interaction of oxygen is shown to strongly contribute 
to the intra-dimer isotropic exchange. The results are in good agreement with those derived from 
experimental data, both regarding the orientation of the Dzyaloshinskii-Moriya vectors and the 
magnitude of all exchange integrals. The microscopic analysis is confirmed by the results of Green 
function's and total energies difference methods. 

PACS numbers: 71.70.Gm, 75.30.Et 



I. INTRODUCTION 

The quasi- two-dimensional compound SrCu2 (603)2 
has attracted a lot of attention over the past ten years 
due to its spin-gapped behavior, its finite magnetization 
below the critical magnetic field deduced from the spin 
gap, and more importantly its magnetization plateaus 
at 1/8, 1/4 and 1/3 of the saturated magnetizationii^i^ 
Based on the quasi- two-dimensional structure of the com- 
pound and on a number of experimental facts, the Hamil- 
tonian expected to provide an accurate description of the 
magnetic properties of that compound must include both 
Heisenberg symmetric exchange processes and antisym- 
metric Dzyaloshinskii-Moriya (DM) interactions. It is 
defined by: 

n.n n.n.n 

+DY,[SixSj]+D' Y,[S,xS,], (1) 

n.n n.n.n 

where J {D) and J' {D') define intra- and inter- 
dimer symmetric (Dzyaloshinskii-Moriya) exchange in- 
teractions. The Heisenberg model obtained when D and 
are set to zero is known as the Shastry-Sutherland 
model^. The susceptibility and the main features of the 
magnetization curve have been interpreted in the context 
of this Shastry-Sutherland model, leading to estimates of 
J and J', while the dispersion of triplet excitations, the 
low-field uniform and staggered magnetizations and ESR 
measurements have provided estimates of the DM inter- 
actions. A very useful complementary source of informa- 
tion is usually given by first-principles calculations, which 
have provided valuable results even in low-dimensional 
quantum spin systems with complex structures. Sur- 
prisingly enough, no ab-initio results have been reported 
so far for SrCu2 (303)2- 



In this paper, we report on the first ab-initio inves- 
tigation of the electronic structure and magnetic prop- 
erties of SrCu2 (603)2- Our results show that an ac- 
curate treatment of the oxygen magnetization is crucial 
for the description of isotropic exchange interactions in 
SrCu2 (603)2- In order to calculate the Dzyaloshinskii- 
Moriya interactions between magnetic moments, we have 
used the microscopic expression derived by Moriya.^ For 
that purpose, the hopping integrals between different 
Wannier orbitals centered at Cu atoms have been ob- 
tained using a projection procedure. 

The paper is organized as follows: in the next section, 
we shortly describe the crystal structure of SrCu2 (603)2 
and present the results of LDA calculations. Sec. HI and 
IV contain the analysis of the isotropic and anisotropic 
exchange interactions. The results of LSDA+U calcula- 
tions are presented in Sec.V. We discuss and briefly sum- 
marize our results in Sec. VI and VII . 



II. RESULTS OF LDA CALCULATIONS 

The simplified crystal structure of SrCu2 (603)2 is pre- 
sented in Fig.l. Each copper atom has one nearest- 
neighbor Cu atom and four next-nearest-neighbor Cu 
atoms in CU6O3 layers. In Sr layers, the distance be- 
tween Cu^+ ions is shorter than in CU6O3 layers. How- 
ever, our results demonstrate, as we will show below, 
that the magnetic interactions inside the CU6O3 layers 
are much stronger than in the Sr layers. 

The electronic structure calculation of SrCu2 (603)2 
has been performed using the Tight 6inding Linear- 
Muffin-Tin- Orbital Atomic Sphere Approximation (T6- 
LMTO-ASA) method in terms of the conventional local- 
density approximation.^ We used the known crystal 
structure data.^ The radii of atomic spheres have been 
set to r(Cu) = 2.2 a.u., r(0) = 1.65 a.u., r(6) = 1.3 a.u.. 
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tion basis. ^^'^^ We have defined the Wannier functions 
centered on the copper orbitals with x^ — y^ symmetry us- 
ing a projection procedure. The resulting orbitals (pre- 
sented in Fig. 5) have a strong contribution from the oxy- 
gen atomic wave orbitals and can be expressed through 
the following linear combination: 



(2) 



where a and P are the amplitudes of the copper and oxy- 
gen atomic wave functions in the Wannier orbital Wi (x) . 



III. ISOTROPIC EXCHANGE INTERACTIONS 



FIG. 1: Simplified three-dimensional structure of 
SrCu2 (603)2. The cycles are copper atoms. The dashed 
lines correspond to inter-dimer interaction paths. 



r(Sr) = 3.5 a.u. and nine types of empty spheres were 
added. 

Fig. 2 gives the total and partial density of states pro- 
jected on the constitutional atoms of the SrCu2 (603)2- 
The following characteristics can be seen from Fig. 2. (1) 
The bands lower than - 15 eV mostly consist of 2s orbitals 
of the oxygen atoms of the (BOa)^" complex. (2) The 
density of states near Fermi level comes mainly from Cu 
and O orbitals. The band structure of SrCu2 (603)2 near 
the Fermi level obtained by LDA calculation is presented 
in Fig. 3. There are four well separated bands. The anal- 
ysis of the partial density of states (Fig. 4) shows that 
the main contribution to these bands comes from copper 
orbitals with x^-y^ symmetry. There is however a contri- 
bution of oxygen 2p states due to the strong hybridization 
of Cu 3d and O 2p states (see Fig. 2 and Fig. 4). In this 
situation, the most natural and simplest way to describe 
the magnetism and to take into account the hybridization 
between copper and oxygen is work in a Wannier func- 



The results of the previous section demonstrate that 
within LDA the SrCu2 (603)2 has a metallic ground 
state. This is a standard LDA problem for transition 
metal oxides and in order to overcome it one should take 
into account on-site correlations in the electronic density 
functional. Despite this problem, a microscopic descrip- 
tion of magnetic couplings is still possible on the basis 
of LDA. In order to demonstrate this, we construct a 
model Hamiltonian describing SrCu2 (603)2 in the Wan- 
nier function basis. This Hamiltonian is set to repro- 
duce only four bands of the full Hamiltonian near the 
Fermi level (Fig. 3). We express the LDA structure re- 
sults in terms of a low-energy, few orbital model Hamilto- 
nian using a projection procedure. Let us first consider 
the one-orbital tight-binding Hamiltonian in the Wan- 



nier functions basis Htb 



^ al^aja, where 



t^- ~^ is the hopping integral between Wannier functions 
centered on the 3da.2_^2 orbitals of ith and jth sites. To 
simplify the analysis we divide the calculated hopping 
integrals into two subgroups: intraplane (Fig. 6) and in- 
terplane (Fig. 7) couplings. The intraplane transfers are 
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FIG. 2: Total and partial density of states of SrCu2 (603)2 
obtained in LDA calculations. 



FIG. 3: Band structure of SrCu2 (603)2 near Fermi level (0 
eV). 
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FIG. 4: Partial density of 3d copper states obtained by LDA 
calculations. The blue dashed, black and green solid lines are 
3d states of copper, 3da,2_^2 orbital and 2p states of oxygen, 
respectively. The zero energy corresponds to the Fermi level. 

^in^ = = 12 meV. In turn the interplane hop- 

2_ 2 2_ 2 

ping integrals have the values: t^^ ^ =^ib ^ =8 meV 

and tfc'^' = t^^"^'= tffy'= if-y' = 12 meV. 

Using these hopping parameters, one can estimate the 
magnetic couplings of SrCu2 (603)2- Since the Wannier 
orbitals have a strong contribution from the wave func- 
tions of the oxygen atoms, it is not enough to use the 

simple formula —ff- . There is an additional ferromagnetic 
contribution originating from Hund's rule intra-atomic 
exchange interaction of the oxygen^, which leads to the 
formula: 

where Nqx is the number of oxygen atoms between cop- 
pers and P is the contribution of atomic wave func- 
tions of oxygen to the Wannier orbitals. IJd and 
are on-site Coulomb interactions of copper atom and 
intra-atomic exchange interaction of oxygen, respectively. 
The former interaction can be estimated within con- 
strained LDA, which gives 11(^=8.4 eV. According to spin- 
polarized LSDA+U calculations^ =1.6 eV. The value 




FIG. 5: Wannier orbitals centered at —y^ orbitals of copper 
atoms which belong to the dimer. 




FIG. 6: Representation of couplings between Atom 1 and 
its neighbours which belong to the plane. The interactions 
between Atom 1 and the light grey sites without numbers are 
negligibly small. 



of the copper magnetic moment obtained by LSDA+U 
calculations is 0.72 /i^. One can estimate through the 
magnetic moment of the oxygen atom in the ferromag- 
netic configuration simulated in LSDA+U calculations 
(Table III): = M(O)/2=0.05. Using Eq.(3) and the 
parameters defined above leads to: J12 = 10.2 meV. The 
value of J12 is larger than that derived from experimental 
data. The problem might be that Eq.(3) is sensitive to 
the form of the Wannier functions. For instance, if 
was equal to 0.054 (instead of 0.05), then J12 would be 
equal to 7.5 meV, in excellent agreement with the exper- 
imental value. A detailed comparison of our results for 
the magnetic couplings with those extracted from exper- 
iments is presented in section VI. 

On the other hand, there is no overlap at oxygen 
atom between Wannier functions which belong to dif- 
ferent dimers and therefore the inter-dimer interactions 




FIG. 7: Interplane interactions between copper atoms. The 
interactions between Atom 1 and the light grey sites without 
letters are negligibly small. 
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can be calculated using the standard expression: 



(4) 



and 



// 



C (^)^p. (^)^l (^O^p. (^0 



dxdx' . 



\x — X' 



(9) 



This leads to Ji3=3.6 meV, Ji7=0.2 meV, Ji9=0.13 meV, 
JiA 0.06 meV and Jic 0.13 meV. 



IV. ANISOTROPIC EXCHANGE 
INTERACTIONS 

In a pioneer investigation, Moriya^ pointed out two 
main contributions to the anisotropic exchange interac- 
tion. The first one is the kinetic Dzyaloshinskii-Moriya 
interaction, which can be expressed in the following form: 



^3 



kin 
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One can estimate these integrals through complex spher- 
ical harmonics and Slater integrals. We found that, due 
to symmetry, both integrals of Eq.(8) and Eq.(9) are 
identically zero. Therefore, we only consider the kinetic 
Dzyaloshinskii-Moriya interaction in the following. 

To perform the microscopic analysis and calculate the 
kinetic Dzyaloshinskii-Moriya interactions, we define the 
hopping parameters of the following general tight-binding 
Hamiltonian that includes five Wannier orbitals centered 
at Cu sites: 



rr5 orb 



(10) 



where n and n' denote the ground state Wannier orbitals 
centered at the ith and jth sites, while t^p and C'^p are 
transfer integrals without and with spin-orbit coupling. 
In the case of SrCu2(B03)2, the ground state Wannier 
function is centered on the 3d2;2_^2 orbital of copper. If 
the energy difference between the ground state and the 
excited states is larger than the spin-orbit coupling, then 
we can treat the spin-orbit interaction as a perturbation. 
In this case C^^ ^ is given by 
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where A is the spin-orbit coupling constant, is the 
matrix element of the orbital angular momentum be- 
tween the mth excited state and the nth ground state 
Wannier functions which are centered at ith ion, while 
represents the energy of the nth Wannier orbital at 
the ith ion. 

The second contribution to anisotropic exchange in- 
teractions is of Coulomb interaction origin,^ and is given 

by 



where /c, = xy^ yz^ 3z^ — r^, xz, x'^ —y'^. Using the pro- 
jection procedure^^ we have calculated the hopping inte- 
grals between the ground state Wannier orbitals of x'^ — y^ 
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The hopping integrals between the ground state {x^ — y'^) 
and excited {xy, yz^ 3z^ — r^, xz) states Wannier orbitals 
are given 
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a kind of inter-site Coulomb exchange interaction which 
is non-diagonal with respect to the orbitals. Assum- 
ing that the Wannier functions of x'^ — y'^ symmetry are 
ground state orbitals, these intersite Coulomb exchange 
interaction integrals can be expressed in the atomic wave 
function basis in the following form: 



// 
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In contrast to the one-orbital model, the value of the 
intra-dimer hopping integral between Wannier orbitals of 
x'^ — y'^ symmetry becomes smaller. This is due to the fact 
that the hybridization between filled and vacant orbitals 
is explicitly taken into account in the five-orbital model. 
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We can estimate the isotropic exchange interactions in 
the case of the five-orbital model using Eq.(3), J12 = 7.8 
meV, Ji3=3.5 meV, J17 = 0.23 meV, J19 = 0.15 meV, 
JiA = 0.13 meV and Jic = 0.09 meV. These values are 
in better agreement with the experimental estimates. 

The energies of the Wannier orbitals obtained by the 
projection procedure^^ are shown in Table I. Using these 



TABLE I: The energies, Cm of the Wannier orbitals obtained 
using the projection procedure (in eV). The zero energy cor- 
responds to the Wannier orbitals of 3z^ — symmetry 

xy yz Sz —r zx x — y 



0.26 0.20 0.20 2.54 

hopping integrals, we find that Eq.(5) can be reduced 
to the following expression for the x and y components 
of the intradimer Dzyaloshinskii-Moriya interaction (i=l 
and j=2): 



r)X 



16A 



yz ^x^ 
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TABLE IL Calculated anisotropic exchange interaction vec- 



tors Dij (in meV). 



(iJ) 



(1.2) (0.25;-0.25;0.0) 

(1.3) (0.0; -0.07 ; 0.1) 

(1.4) (-0.07; 0.0; -0.1) 

(1.5) (0.0; 0.07; 0.1) 

(1.6) (0.07;0.0;-0.1) 



and 



D and correspond to D12 and Dij (where j =3,4,5,6), 
respectively. Based on Eq. (11-12), we can conclude that 
the source of intra-dimer anisotropic exchange interac- 
tion is the hopping process between Wannier orbitals of 
x'^ — y'^ and yz{xz) symmetry. The microscopic origin of 
the inter-dimer Dzyaloshinskii-Moriya interaction is the 
transfer process between x'^—y'^ and xy Wannier orbitals. 
It is interesting to note that there are x and y compo- 
nents of inter-dimer anisotropic couplings. The compari- 
son with experimental data is presented in Sec.VL In the 
next section we present the results of LSDA+U calcula- 
tions, which give additional confirmation of our model 
consideration for the isotropic interactions. 



DI2 = -iTTT^ -tf,-y\if-y\ (12) 



RESULTS OF LSDA-hU CALCULATIONS 



The z component of the inter-dimer (i=l and j=3,4,5,6) 
coupling is given by 
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We are now in a position to calculate Dzyaloshinskii- 
Moriya interactions using Eq.(ll-17) and the value of 
spin-orbit coupling A = 0.1 eV. These results are pre- 
sented in Table II and Fig. 8. In the notation of Eq.(l) 



The results of LDA calculations have demonstrated 
that the magnetic properties of SrCu2 (603)2 can be 
reproduced correctly within the microscopic model ap- 
proach based on the LDA calculations. However, LDA 
cannot correctly describe the insulating ground state of 
SrCu2 (303)2- In order to overcome this problem, we 
have used the LSDA+U approximation^^ which takes 
into account the Coulomb correlations of localized states 
neglected in LDA. The effective Coulomb interaction U^^ 
and the effective intra-atomic exchange of the cop- 
per atoms, which represent external parameters in a self- 
consistent cycle of the LSDA+U scheme, are determined 
from the first-principle calculation by constrained LDA. 




FIG. 8: Schematic representation of symmetry of 
Dzyaloshinskii-Moriya interactions. 



TABLE III: Results of LSDA+U calculations for ferromag- 
netic and antiferromagnetic configurations. E^^p is the value 
of the energy gap (in eV). M(Cu) and M(0) are the values 
of the magnetic moments of copper and oxygen atoms which 
are located between the copper atoms in dimer (in /hb)- Etot 
is the relative total energy of the system in different magnetic 
configurations (in meV). 



M(Cu) M(0) Et 



FM 
AFM 



2.2 
2.2 



0.72 
0.72 



0.1 
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The calculation scheme has been described elsewhere. 
The Coulomb interaction parameter Ud and the intra- 
atomic exchange have been estimated to be Ud=9A 
eV and Jf =1 eV. 

We have performed LSDA+U calculations for two mag- 
netic configurations which are presented in Fig. 9. The 
results of these calculations are presented in Table IIL In 
both cases, SrCu2(B03)2 is an insulator with an energy 
gap of 2.2 eV. The antiferromagnetic configuration has a 
lower energy. 

Let us describe the inter-dimer interactions. For that 
we use the Green's function method. Following Lichten- 
stein et al}^ ^ we determine the exchange interaction pa- 
rameter between copper atoms via the second variation 
of the total energy with respect to small deviations of the 
magnetic moments from the collinear magnetic configu- 
ration. The exchange interaction parameters J^j can be 
written in the following form:— 



Jii = - delm (A—' G^'f" Af'-'" 

J —CO / 

m" ^m'" 

where m is the magnetic quantum number, the on-site 
potential A^^^' = H^"^' - H^f and the Green's func- 
tion is calculated in the following way 



(18) 



is a component of the nth eigenstate, 
and is the corresponding eigenvalue. This 

leads to the following inter-dimer exchange interac- 
tions: Ji3=Ji4=-^i5 =-^16=4.1 meV, Ji7=Ji8=0.17 meV, 
-^i9=-^iio = Jill = Jii2 0.06 meV, Jia=Jib 0.06 
meV, Jic = JiD = JiE=JiF= 0.02 meV. 

One can see that in the case of the antiferromagnetic 
configuration there is no magnetic moment at the oxygen 
atom (Table III). But in the ferromagnetic configuration, 
the oxygen has a small moment. This fact supports the 
scenario proposed in Ref. [5]. Therefore, to calculate the 
intra-dimer exchange interaction, it is necessary to take 
into account the change of oxygen magnetization. One 
can do this using the method of total energies difference. 
For that purpose, we write the Heisenberg Hamiltonian 




FIG. 9: Ferromagnetic and antiferromagnetic configurations 
simulated in LSDA+U calculations. 



in the following form: 

H — ^ JijSiSj. 



(19) 



1>J 



The total energies of the ferromagnetic and antiferromag- 
netic configurations presented in Fig. 9 are 



EpM = J12S -\-4JisS -\-2JiyS 



and 



(20) 



(21) 



Eafm = —J12S — 4J13S —2JnS . 
Therefore, the dimer exchange interaction J12 is given by 
EpM — Eafm — SJisS^ — iJnS^ 



J12 = 



252 



(22) 



Using calculated values of the inter-dimer couplings ( J13 
= 4.1 meV and J17 = 0.17 meV) and the values of the 
total energies from Table III, one can calculate the dimer 
exchange interaction for S=l/2, which leads to J12 = 7.3 
meV. 
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VI. COMPARISON WITH EXPERIMENT 

The exchange interactions obtained with different 
methods are summarized in Table IV, together with the 
values most often used to explain experimental data. The 
Green's function method is in principle the most accu- 
rate, and if its results are in good agreement with those 
obtained by LDA using microscopic exchange formulae, 
which is the case here, the results are expected to be quite 
accurate. The agreement with experimental estimates is 
quite good for J12 and J13. The inter-plane coupling 
JiA is however predicted to be about 10 times smaller 
than what has been suggested from fitting the suscepti- 
bility, a conclusion to be taken seriously in view of the 
overall good agreement regarding the other parameters. 
Besides, the dominant couplings beyond J12 and J13 are 
predicted to be the in-plane exchange couplings J17 and 
Ji9, parameters which have not been considered so far in 
theoretical models of SrCu2(B03)2. 

Let us go beyond these general statements and look 
in more details at how well the calculated exchange pa- 
rameters can reproduce using basic experimental facts 
such as the spin gap and the temperature dependences 
of the magnetic susceptibility and specific heat. Let us 
start with the spin gap analysis. According to different 
experiment^, the value of the spin gap varies from 2.6 
meV (nuclear quadrupole resonance) to 3.1 meV (nuclear 
magnetic resonance). Based on the calculated exchange 
interactions (Table IV), if we estimate the spin gap using 
the following expression^ 

^ = w-{^f-\{'-^f-\{^n (23) 

J12 ^ J12 o J12 

the spin gap value is 9.1 meV in the case of the one- 
orbital model. This is about three times larger than ex- 
perimental value. A better agreement is achieved in the 
five-orbital model and LSDA+U calculations, A^''''^=5.8 
meV and A^'^^^+^=4.2 meV. Let us note however that 
J\zl J\2 — 0.6 is not that small, and higher order correc- 
tions to the estimate of Eq. ([23|) are expected to further 
reduce the gap, hence to improve the agreement with the 
experimental determination of the gap. 

The results of measurements of magnetic susceptibility 
and specific heat can be reproduced using the Heisenberg 
model with different sets of isotropic parameters. For 
instance, in the review of Miyahara and Uedai, there 
are five different sets of exchange couplings which vary 
from J — 6.6 meV and J' = 4.1 meV to J = 7.3 meV 
and J' — 4.6 meV. The latter J and J' correspond to the 
best fitting to the experimental data. These values are in 
reasonable agreement with our exchange interactions for 
the five-orbital model (J12 = 7.8 meV and Ji3= 3.5 meV) 
and in good agreement with those obtained by LSDA+U 
calculations (J12 = 7.3 meV and Ji3= 4.1 meV) 

The presence of a finite magnetization well below the 
expected critical field for the gap closing is a manifes- 
tation of the fact that triplet states are mixed into the 
ground state. This is an effect of the Dzyaloshinskii- 



TABLE IV: Values of exchange interactions J^j between mag- 
netic moments of SrCu2(B03)2 system (in meV). 

LDA^ LDA' LDA+U^ Exp, (after Ref. \ ) 



J12 


10.6 


7.8 


7.3 


Jl3 


3.6 


3.5 


4.1 


Jl7 


0.2 


0.23 


0.17 


Jl9 


0.13 


0.15 


0.06 


JlA 


0.06 


0.13 


0.06 


Jic 


0.13 


0.09 


0.02 



" One-orbital model. 
' Five-orbital model. 
^ Green's function approach. 

Moriya interaction. The directions (Fig. 8) and values 
(Table II) of calculated Dzyaloshinskii-Moriya vectors 
agree well with those obtained from the interpretation 
of neutron, ESR and NMR experiments. ^^^^ For instance, 
the intra-dimer anisotropic vector is perpendicular to the 
bond direction. The inter-dimer Dzyloshinskii-Moriya in- 
teractions lie mainly along the z axis. The calculated 
ratios J = 0.032 and \D'\/J = 0.016 (with J = 7.8 
meV) are in good agreement with those deduced from 
experiments {\D\/J = 0.034 and \&\/J = 0.02). 

VII. CONCLUSION 

In this paper we have presented the results of an 
ab-initio investigation of the magnetic properties of 
SrCu2(B03)2. It is found that the ferromagnetic contri- 
bution of the intra-atomic exchange interactions of oxy- 
gen plays a crucial role to account for the intra-dimer 
isotropic exchange interaction. The microscopic origin of 
Dzyaloshinskii-Moriya interactions has been analyzed us- 
ing calculated hopping integrals between different Wan- 
nier orbitals of copper atoms. In general, there are two 
contributions to anisotropic exchange interaction which 
are kinetic and Coulomb Dzyaloshinskii-Moriya interac- 
tions. In the case of SrCu2(B03)2 the Coulomb contri- 
bution is zero by symmetry. However, we expect that 
the latter plays important role in the case of more dis- 
torted frustrated cuprates such as Cu2Te2 05(Br,Cl)2 and 
(Tl,K)CuCl3. 
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